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Abstract 

New limits on the Higgs mass from LEP and new calculations on the radiative (pen- 
guin) decay of the 6-quark into 57 restrict the parameter space of the Constrained Minimal 
Supersymmetric Standard Model (CMSSM). 

We find that for the low tan [3 scenario only one sign of the Higgs mixing parameter is 
allowed, while the high tan (3 scenario is practically excluded, if one requires all sparticles 
to be below 1 TeV and imposes radiative electroweak symmetry breaking as well as gauge 
and Yukawa coupling unification. For squarks between 1 and 2 TeV high tan (3 scenarios 
are allowed. We consider especially a new high tan (3 = 64 scenario with triple unification 
of all Yukawa couplings of the third generation, which shows an infrared fixed point 
behaviour. 

The upper limit on the mass of the lightest Higgs in the low (high) tan (3 scenarios is 
97±6 (120±2) GeV, where the largest error originates predominantly from the uncertainty 
in the top mass. 



1 Introduction 



Recently several new results have appeared, which have a considerable impact on the allowed 
parameter space of the Constrained Minimal Supersymmetric Model (CMSSM). As constraints 
we consider gauge unification and radiative electroweak symmetry breaking assuming the su- 
pergravity inspired scenario of the soft breaking of the SUSY masses with common scalar and 
gaugino masses at the GUT scale The additional reduction of the parameter space by relic 
density, b — r Yukawa coupling unification and the new Higgs limits from LEP as well as the 
present b — » sj rate is calculated. 

The most significant new results are: 

• the SM Higgs mass limit is now 89.3 GeV at 95% C.L. for the SM Higgs [§, which 
corresponds to the CMSSM case too, since all the other Higgses are heavy and decouple 
in the CMSSM limit. 

• Higher order corrections to the b — > calculation became available P, |], [|] . As exper- 
imental value we use the combined value from the CLEO Coll.|J and ALEPH Coll.j7| 
B1Z{b — > 57) =(2.54 ± 0.56) • 1CT 4 , which is consistent with the SM expectation: BlZ{b — > 
57) =(3.6 ± 0.3) ■ 10~ 4 . This SM value was calculated for a s = 0.122, which is the best 
value from the Z° data of the four LEP experiments^. Since the calculations are usually 
quoted as the ratio of the b — > 57 rate and the semileptonic decay rate of the b-quark, one 
needs for absolute predictions BlZ(b — > X c en), for which we use 10.5 ± 0.4% 0. 

In a recent paper the theoretical uncertainties on the photon spectrum are studied in 
next-to-leading order as well0. They could increase the CLEO value by as much as 14%. 
Such a correction improves the agreement between data and the SM prediction. Given 
the present uncertainty of 22% in the published data, this correction has not been taken 
into account, so we stick to the published values and added in the analysis the theoretical 
uncertainty from the SM prediction, the semileptonic branching ratio uncertainty and the 
experimental errors in quadrature. 

Our statistical analysis, in which the probability for each point of parameter space is deter- 
mined by a x 2 fit, has been described previously ]TD[. The masses of the third generation and the 
couplings constants are taken from LEP results presented at the Jerusalem Conf . || , except for 
the top mass from the combined data from DO and CDF, which was updated: mt = 173.9 ±5.2 
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Since the common GUT masses m and mi/ 2 for the spin and spin 1/2 particles are 
strongly correlated, we perform the fit for all combinations of mo and mi/2 between 100 GeV 
and 1 TeV in steps of 100 GeV. 



2 Results 

2.1 Constraints from the top and bottom mass 

The requirement of bottom-tau Yukawa coupling unification strongly restricts the possible 
solutions in the m t versus tan/3 plane, as discussed in our previous paper[|TU . 

Only three values of tan/5 give an acceptable x 2 fit for m t = 174 GeV, as shown in the 
bottom part of fig. [T]. Clearly, the value tan f3 = 64 gives a worse \ 2 than the values at 35 and 
1.65, which is mainly due to the fact that this solution cannot fit simultaneously the top and 
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bottom mass with a single tan [3 if b — r unification is required ( either m t or m^ too large by 
about 3a). 



Therefore it was not considered in our previous paper |L0|. We consider it here, since the 
solution with tan f3 ~ 35 does not give a better x 2 anymore in the combined fit of all constraints 
because of the smaller errors from the gauge couplings, mt and b — > sj. 

Furthermore, this solution has the unique feature of a possible triple Yukawa unification: 
all 3 Yukawa couplings are driven to an approximate fixed point, which yield approximately 
the correct masses of the leptons and quarks of the third generation. Since the mass deviation 
are at the 3a level, we consider it here, especially since the solution with tan (3 ~ 35 does not 
give a better \ 2 with the smaller errors from the gauge couplings, mt and b — > 37. The fixed 
point is reached for a large range of initial values of the Yukawa couplings at the GUT scale, 
as shown in fig. |2|. 

The difference between the two solutions at high tan/3 corresponding to opposite signs of 
/j, stems from finite loop corrections to the bottom quark mass involving squark-gluino and 
stop-chargino loops. These corrections are small for low tan [3 solutions, but can become as 



large as 10-20% for the high tan/3 values[T(|, since the dominant corrections are proportional 
to /itan/3. Consequently they change sign for different signs of /1. As shown in the middle 
part of fig. [I] the Yukawa couplings become of the same order of magnitude at large tan (3. For 
tan (3 = 64 Yj, varies so rapidly that it is easy to find a solution where all Yukawa couplings are 
equal, but large at the GUT scale. The low energy value is then approximately given by a fixed 
point solution for each of the 3 Yukawa couplings, i.e. the solutions of the RGE are insensitive 
to the starting values at the GUT scale. 

Radiative electroweak symmetry breaking is only possible in that case, if one assumes non- 
universal values for the mass parameters in the Higgs potential. 

As will be discussed below, the b — > 57 rate is also sensitive to the sign of \x and it turns out 
to be difficult to have solutions, where both b — r unification and the b — » 57 rate are correctly 
described. 



2.2 b — > S7 decay rate 

Recently the Next-to-Leading Order (NLO) QCD contributions to the b — > sj decay rate have 
been calculated in the SM||. They increase the SM value by about 10%. The NLO electroweak 
(EW) contributions decrease the result by a few [3J]. In NLO the large theoretical uncer- 
tainty from the scale dependence is now reduced to the 10% level (see e.g. ref. || for a recent 
summary on the discussion of the uncertainties). The NLO corrections can be simulated in the 
lowest order calculation by choosing a renormalization scale fi ~ 0.6m&. Since the NLO have 
not yet been calculated for the MSSM, we will stick to the first order calculations, but choose 
/i = 0.6mf,. Since the NLO corrections are small compared with the experimental uncertainty 
of 22%, we expect that with the present large experimental uncertainties the NLO corrections 
in the MSSM will not change the conclusions significantly. 

The b — > 57 branching ratio is very sensitive to the presence of SUSY particles. Their effect 
has been calculated by Bertolini et al [O where a complete but cumbersome matrix notation 



for the amplitudes was used. We use the more explicit expressions by Barbieri and Giudice [fl3 
However, in the latter paper the intergenerational mixing in squark mass matrices is omitted, 
which is not justified in all regions of parameter space: we find it can change the charged 
current amplitude with charginos in the loop by ~ 10% or more. So the usual assumption that 
intergenerational mixing predominantly introduces only the usually negligible flavour changing 
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Figure 1: The top quark mass as function of tan/5 (top) for values of mo, mi/2 ~ 1 TeV. The 
curve is hardly changed for lower SUSY masses. The middle part shows the corresponding 
values of the Yukawa couplings at the GUT scale and the lower part the x 2 values. If the top 
constraint (m t = 174 ± 5, horizontal band) is not applied, all values of tan/3 between 1.2 and 
70 are allowed, but if the top mass is constrained to the experimental value, only the regions 
tan/3 = 1.65 ± 0.3, tan/3 ~ 35, and tan/3 ~ 64 are allowed. 

neutral current with gluinos and down- type quarks in the loop, is not valid. Therefore we 
extended the explicit formulae of Barbieri and Giudice with the intergenerational mixing, thus 
obtaining an efficient code, which is important for numerical x 2 analysis. 

At low tan/3 the b — > rate is close to its SM value for most of the plane (see fig. |3]). The 
charginos and/or the charged Higgses are only light enough at small values of m and m 1( / 2 to 
contribute significantly. 

However, for large tan/3 the b — > s'j rate strongly depends on the particle spectrum and 
can be larger than the SM value by factors 3 or more, as shown on the right hand side of fig. 
|3|. The chargino contribution, which is proportional to /i tan /3, changes sign for different signs 
of (1. Through interference with the large operator mixing term it can lower or increase the 
total b — > sj rate. The gluino contribution is noticeable only if mi/2 is small. However, in this 
region the charginos are light too, so the chargino contribution is always larger by an order of 
magnitude. 

For negative [i values the b — > sj rate is usually too large, except if all sparticles are very 
heavy and don't contribute (see fig. right hand side). For positive \i one can reach the values 
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Figure 2: The running of the Yukawa couplings in case Y t = Y& = Y T at the GUT scale (S'O(IO) 
type solution). One can clearly see the approach to the three different fixed points, i.e. the 
low energy value is largely independent of the GUT scale value. Consequently the GUT scale 
values can be choosen to be equal (triple unification). The fixed point values at low energy 
yield correct masses for bottom and tau for tan (3 ps 64; the fixed point of the top mass yields 
= (47r) 2 Y t v 2 sin 2 /3=184 GeV, which is about 2a above the experimental value. The values 
of the Yukawa couplings at the GUT scale are shown in the maximum allowed range: higher 
values of Y t tend to produce negative Higgs masses, while lower values yield a too low top mass 



of b — > sj rate close to the experimental value in most of the parameter region, as is apparent 
from the x 2 distribution in the left hand bottom side of fig. However it is difficult to fit 
simultaneously m&, if the Yukawa coupling YJ, is determined from m T via b — r unification (see 
left hand top part of fig. |J). Note that for this sign of fi the mass is not just slightly wrong, 
for which threshold corrections could be blamed, but for our case off by more than 7a, where 
we used the conservative Particle Data Book errors for the running mass of mj, = 4.2 ± 0.15 
GeV. Recent determinations from Ypsilon mass spectra quote a 50% lower error [14|. 

On the right hand side of fig. |] the sign of fi is reversed. In this case mb can be fitted 
very well, but b — > 57 cannot be fitted for mi/2 < 600 GeV. The reason is that in this case the 
negative chargino loop contribution oc /itan (3 yields a strong positive contribution through the 
interference with the negative operator mixing term, so BlZ(b — > sj) comes out too large, if the 
charginos are light. On the contrary, the large mj, contributions from g — q and — t loops 
proportional to /itan (3 now lower m& sufficiently 

Our conclusions may seem to contradict those of refs. [15| and [16 



In ref. [15] the discrep- 



ancy between the left bottom and left top of fig. |] is solved by minimizing the corrections to 
rrib proportional to a s /xtan/3, which can be done by a small value of \i and a small value of a s . 
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tan (3 = 1.65, (X > 



tan P = 35, |i< 




Figure 3: The predicted b — > 57 branching ratio for low and high tan (5 in units of 10 4 as function 
of m and mi/ 2 . At low tan/3 the prediction is close to the SM value (3.5) in the whole plane, 
while at high tan /3 it can be enhanced by a large factor for light sparticles (mainly charginos). 



Giving up exact gauge unification by using larger errors on the gauge couplings and allowing 
for an offset of a few % for the strong coupling a s at the GUT scale, they find acceptable fits. 

Our b — > sj calculations are in reasonable agreement with the results from ref. [HJ for both 
signs of [x. However, since they do not impose b — r unification, their preferred solution is \x > 0, 
which corresponds to the left bottom corner of fig. f| and ignoring the left top corner. 



2.3 Excluded parameter space and the mass spectrum 

In fig. |^ the total x 2 distribution is shown as a function of m and mm for the three values 
of tan/3 determined above from b — r unification. The areas at low mO and high rai/2 are 
excluded by the LSP constraint, since in that case the lightest f can become the LSP. If Im- 
parity is conserved, a charged LSP is not allowed, since the vacuum would be filled with charged 
relics from the Big Bang. 

The relic density constraint excludes m > 350 GeV in case of small tan/3, as discussed 
previously For large tan/3 the Higgsino mixture of the LSP allows a fast enough decay via 
s-channel Z° exchange, so in that case the requirement Qh 2 < 1 is easily fulfilled. 

As discussed in the previous section the combined requirements of the correct b — > 37 rate 
and b — r unification exclude quite a region of parameter space for large tan /3, as shown by the 
contours in the lower part. 

One observes X 2 minima at m ,m 1/2 around (200,500), (1000,900), and (500,500) for the 
different tan (3 values, respectively, as indicated by the stars. 

Note that the squarks and gluinos are typically above 1 TeV for the high tan/3 solutions. 
Furthermore, the minimal x 2 values are not excellent for high tan/3 : for tan/3 = 64 x^ in =6.1 
from the fitted top mass (m t = 189 GeV), while for tan/3 = 35 x^ n =4.3 from b — > 57 . 
All other x 2 contributions are negligible. For tan/3 = 1.65 Xmin = basically from the 
BlZ(b — ► 57) constraint alone. 
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tan p = 35, jll > 



tan p = 35, jll < 




Figure 4: The \ 2 contributions from the bottom mass and b — > for different signs of /x in the 
high tan/3 solution. Note that b — > 57 is perfect for /j, < 0, while the bottom mass is perfect 
for fj, > 0. Requiring both yields a low x 2 only for heavy sparticles (mi/2 > 650 GeV). 



Apart from the heavy spectra for large tan/3 , one has the problem that the Born term 
Higgs masses are strongly negative, as can be anticipated from the fast running of the soft 
mass terms of the two Higgs doublets, m\ and m 2 ., which receive negative radiative corrections 
proportional to the Yukawa couplings. For large tan/3 both Y t and are large (see fig. |T|). 
Then the pseudoscalar Higgs mass, (m 2 + m 2 , at the Born level), and the lightest Higgs get 
negative masses at the Born level, which can only become positive by large radiative corrections. 
In this case it is extremely important to include the full 1-loop radiative corrections, i.e. one 
has to include both the sparticles and particles in the loops. For the full 1-loop corrections we 
used the formulae from ref. [[17]]. For the lightest Higgs the dominant second order corrections 



from ref. []T8| have been included as well. 

For low tan/3 the present Higgs limits severely constrains the parameter space, as can be 
seen from fig. || which shows the excluded regions in the (7710,7711/2) plane for different signs of 
fi. As mentioned in the introduction the SM Higgs limit of 89.3 GeV is valid for the low tan/3 
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Figure 5: The ^-distribution for the low and high tan (3 solutions. The different shades in 
the projections indicate steps of Ax 2 = 4, so basically only the light shaded region is allowed. 
The stars indicate the optimum solution. Contours enclose domains excluded by the particular 
constraints used in the analysis. 



scenario (tan j3 < 4) of the MSSM too. As is apparent from fig. ^| this limit clearly rules out 
the /i < solution, in agreement with other studies|21|]. However, this figure assumes mt = 175 
GeV. The top dependence on the Higgs mass is slightly steeper than linear in this range (see 
fig. 0). Adding about one a to the top mass, i.e. m t = 180 GeV, implies that for the contours 
in fig. ^| one should add 6 GeV to the numbers shown. Even in this case the \x < solution is 
excluded for a large region of parameter space. Only the small allowed region with m 1( / 2 > 700 
GeV is not yet excluded for m t = 180 GeV. Note that in this region the squarks are well above 1 
TeV, so in this case the cancellation of the quadratic divergencies in the Higgs masses, which is 
only perfect if sparticles and particles have the same masses, starts to become worrying again. 

In summary, if all squarks are required to be below 1 TeV the \x < solutions are excluded 
for the low tan (3 scenario, even if the uncertainty from the top mass is taken into account. For 
H > practically the whole plane is allowed, except for the left bottom corner shown on the 
top right hand side of fig. |6|. 
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Figure 6: Contours of the Higgs mass (solid lines) in the m , my 2 plane (above) and the Higgs 
masses (below) for both signs of \i for the low tan/3 solution tan/? = 1.65 for M t = 175 GeV. 

The upper limit for the mass of the lightest Higgs is reached for heavy squarks, but it 
saturates quickly, as is apparent from the bottom row in fig. || For m = 1000, mi/ 2 = 1000, 
which corresponds to squarks masses of about 2 TeVQ, one finds for the upper limit on the 
Higgs mass in the low tan f3 scenario: 



max 
m h 



97 ±6 GeV, 



where the error is dominated by the uncertainty from the top mass. If one requires the squarks 
to be below 1 TeV, these upper limits are reduced by 4 GeV. 

It agrees well with the value from Casas et al.: = 97 ± 2 GeV ||20|| . In both analysis the 
Renormalization Group Equations are used to determine the trilinear coupling at low energies 
and n from EWSB, so the mixing in the stop sector is fixed, once the sign of \i is choosen. 
Furthermore in both papers solutions close to the infrared fixed point are considered, which are 
required in our case from the EWSB. Our error on the upper limit is larger than the one from 
Casas et al., since they did not consider the error on the top mass. 

For high tan/3 the upper limit on the Higgs mass in the CMSSM is: 



m. 



120 ± 2 GeV. 



The error from the top mass is small, since the high tan f3 fits prefer anyway top masses around 
190 GeV. 



1 Explicit analytical expressions for the sparticle masses as function of the SUSY parameters can be found in 
ref. M- 
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Figure 7: The mass of the lightest CP-even Higgs as function of the top mass at Born level 
(dotted lines), including complete one-loop contributions of all particles (dashed lines). Two- 



loop contributions |18[ reduce the one- loop corrections [T7| significantly as shown by the shaded 
area (the upper boundary corresponds to fj, > 0, the lower one to n < 0). The solid line just 
below the dashed line is the one-loop prediction from the third generation only, which apparently 
gives the main contribution. The upper scale indicates the value of tan/3, as calculated from 
the top mass by the requirement of b — r-unification. 

3 Summary 

It is shown that the CMSSM can fit simultaneously the constraints from 

• gauge couplings unification; 

• b — r Yukawa couplings unification; 

• radiative EWSB; 

• life time of the Universe. 

The present Higgs limit of 89.3 GeV largely excludes the fi < solution for the low tan/3 
scenario. The fi > solution is just within reach of the coming LEP II runs above 189 GeV, 
since the upper limit on the Higgs mass is 97 ± 6 GeV. Here the central value of 97 GeV was 
calculated for m t = 175 GeV and = 103 GeV is reached for m t = 180 GeV. If one requires 
the squarks to be below 1 TeV, these upper limits are reduced by 4 GeV. 

The maximum value of the lightest Higgs mass in the high tan /3 scenarios is 120 ± 2 GeV, 
which is beyond the reach of LEP II. However, the high tan j3 solutions have serious finetuning 
problems, because of the fact that all Yukawa couplings are large, which cause the Higgs masses 
at tree level to be strongly negative (typically -1 TeV), so the radiative corrections have to be 
positive and very large to offset this large negative "starting" value. Furthermore, the solutions 
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with minimal \ 2 require squarks above 1 TeV, which have an additional finetuning problem 
because of the non-cancellation of the quadratic divergencies to the Higgs masses. These 
quadratic corrections only cancel perfectly, if sparticles and particles have the same mass. 



Note added in proof. 

After completing this work, first results on the NLO corrections to the b — > 57 rate in the MSSM 
became available ||22||. They change the amplitudes by a similar amount as the SM amplitudes, 



but due to the large cancellations in particlular regions of parameter space, the effect can be 
large sometimes. This will be studied in a future paper. We thank G.F. Giudice for a discussion 
on this point. 
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